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NONEQ U ILilBRi UM  STRUCTURE  OF  HYDROMAGNET IC 
GAS-IONIZING  SHOCK  FRONTS  IN  ARGONt 

by 

Martin  I.  Hoffert* 

Polytechnic  Institute  of  Brooldyn 

SUMMARY 

This  study  deals  analytically  with  the  structure  of  gas -ionising  hydromagnetic 
shock  waves.  Since  these  waves,  by  definition,  must  have  non-electrically-conducting 
upstream  states,  their  existence  at  very  high  shock  temperatures  must  be  ruled  out  on 
the  physical  grounds  that  forward -radiated  precursor  ionization  makes  the  unshocked 
gas  conducting.  A  "low  temperature"  collisionally-ionizing  shock  with  ol’ique  magnetic 
field  is  studied  here  to  determine  whether  certain  concepts  which  exist  in  the  current 
literature  are  relevant.  Nondimensionalized  equations  governing  the  nonequilibrium 
structure  of  such  a  front  propagating  into  un-ionized  argon  are  formulated  using 
ionization  rates  and  an  electron  energy  equation  developed  in  an  earlier  paper. 

Comparison  of  the  magnitudes  of  viscous  and  magnetic  Reynolds  numbers  within  this 
front  indicates  that,  if  a  structure  exists,  it  must  consist  of  a  narrow  "imbedded" 
viscous  shock  standing  upstream  of  a  much  wider  hydromagnetic  interaction  and  oniza- 
tion  relaxation  zone.  Hence,  a  modified  form  of  the  Zeldovich-von  Neumann-Daring 
(ZND)  approximation  is  applic-ble  to  the  structure  problem.  It  is  shown  that  :.n  this 
approximation  nontrivial  steady-state  structures  cannot  be  constructed  for  "fe^t"  gas- 
ionizing  shocks.  On  the  other  hand,  solutions  are  possible  for  "slow"  waves,  and 

\  'v‘“- 

The  research  has  been  conducted  under  Contract  Nonr  839(38)  for  PROJECT  DEFENDER, 
and  was  made  possible  by  the  support  of  the  Advanced  Research  Projects  Agency  under 
Order  No.  529  through  the  Office  of  Naval  Research. 

^Graduate  Student;  also,  Senior  Scientist,  General  Applied  Science  Laborato  ies,  Inc., 
Westbury,  Hew  York. 


these  are  obtained  numerically  for  a  femily  of  hydromagnetically  oblique  shocks  at 
Mach  number  Mi  =  20  and  Alfv&i  number  M^  =  1//Z  with  parametrically  varied 
values  of  the  upstream  electric  field.  In  contrast  to  previous  expectations,  the 
upstream  electric  field  is  not  uniquely  defined  by  the  structure.  Because  the  slow 
solutions  are  effectively  exothermic,  to  the  point  where  their  post-shock  temperatures 
are  associated  with  radiation -induced  precursor  ionization,  it  seems  likely  that  caly 
the  solution  with  the  upstream  electric  field  corresponding  to  a  pure  hydromagnetic 
shock  has  physical  significance. 
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1.  INTRODUCTION 


In  recent  years,  a  number  of  'nvestigators  have  contributed  to  the  formula¬ 
tion  of  a  theoretical  model  descriptive  of  the  so-called  gas-ionizing  hydromagnetic 
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shock  wave  (Xulikovskii  and  Lyubimov  ,  Kunkel  and  Gross  ,  Helliwell  ,  Chu  , 

Woods^,  May  and  Tendys1^  and  Taussig**’  *^).  These  waves  are  thought  to  exist, 
for  example,  in  electromagnetic  shock  tubes.  All  the  c.forementioned  authors  either 
postulate  or  imply  that  the  structure  of  these  waves  conforms  to  the  folio  /ing 
archetype  (this  d  scription  can  also  be  taken  as  a  definition  of  a  "gas -ionizing 
hydromagnetic  shock"  in  the  present  context):  Upstream,  the  gas  is  un-ionized, 
electrically  non-conducting  and  hence  uncoupled  from  the  magnetic  fluids  through 
which  the  shock  mc«res.  Consequently,  the  leading  edge  of  the  front  develops  pre¬ 
cisely  as  an  ordinary  hydrodynamic  shock.  Because  of  collisional  ionizing  reactions 
associated  with  the  rising  temperature,  an  electrically  conducting  (hence  hydro- 
magnetically  active)  plasma  is  created  somewhere  in  the  shock  interior.  It  follows 
that  the  overall  structure  is  hybrid  in  nature,  being  partly  hydrodynamic  and  partly 
hydromagnetic. 

The  most  distinctive  implication  of  this  archetype  is  that  the  Rankine- 
Hugoniot  l  onditions  are  no  longer  sufficient  to  predict  the  downstream  state  of  the 
shock  in  terms  of  the  upstream  state  and  the  shock  velocity.  This  is  because,  in 
contrast  to  purely  hydromagnetic  shocks,  the  upstream  gas-frame  electric  field  is 
not  uniquely  defined  in  terms  of  the  upstream  velocity  and  magnetic  fields:  As  a 
non-electrical  conductor,  the  unshocked  gas  is  incapable  of  sustaining  a  current 
flow,  so  the  upstream  boundary  condition  of  no  currents  in  the  undisturbed  gas  is 
automatically  satisfied  for  any  electric  field. 

It  has  been  argued  that  an  analytic  prediction  of  the  electric  field  requires 
an  analytic  and  physically  correct  solution  for  the  ionizing  wave  structure*^.  In 
order  to  gain  some  insight  into  the  structure  problem,  prior  studies**  ***’  *^  have 
assumed  temperature-dependent,  step-function  models  for  the  variation  of  electrical 
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conductivity  7  within  the  shock,  i.  e. ,  0=0  for  T'<T*'  ando^O  for  T'>T*',  where  T' 
is  the  gas  temperature  and  T*'  is  some  "reference"  temperature.  Analysis  by 
May  and  Tendys1^  indicates  that  shock  structure  integral  curves  deduced  from  this 
mode1  are  applicable  only  when  T*'  is  of  the  same  order-of -magnitude  as  u.e 
characteristic  (first)  ionization  temperature  of  the  unshocked  gas. 

The  present  study  is  concerned  with  obtaining  shock- structure  solutions 
(if  any  exist)  which  arc  consistent  with  the  gas -ionizing  archetype  and  which  also 
incorporate  realistic  representa  tons  of  transport  and  rate  processes  in  a  colJi- 
sionally -ionizing  monatomic  ga  t ,  argon  in  particular.  It  is  motivated  by  a  realiza¬ 
tion  that  by  misrepresenting  the  physicr  of  high  temperature  gases,  the  step-function 
temperature -dependent  conductivity  approach  can  give  qualitatively  misleading  results 
for  two  different  reasons:  (1)  If  the  internal  shock  temperature  approaches  the 
ionization  temperature,  as  May  and  Tendys  suggest,  tho  gas  becomes  fully  ionized 
almost  immediately  (since  each  interparticlc  collision  has  enough  energy  on-the- 
average  to  "knock  off"  ar.  outer  electron),  but  at  these  temperatures  radiation- 
induced  precursor  ionization  levels  are  sufficiently  high  so  as  to  preclude  any 
reasonable  interpretation  of  the  upstream  state  as  un -ionized  .  Consequently,  the 
"gas -ionizing"  archetype  with  its  implied  electric  field  indeterminacy  is  violated 
and  the  shock  is  not  gas-ionizing,  in  the  present  context.  (2)  Another  possibility, 
the  one  actually  explox  sd  in  this  paper,  is  that  of  a  "low  temperature"  gas-ionizirp, 
shock,  i.e.,  a  front  creating  a  nonequilibrium  pli  jma  in  which  photo-ionization  is 
realistically  negligible  compared  to  collisional  ionization.  In  this  latter  case  the 
conceots  embodied  in  the  archetype  may  still  be  revelant  but  the  step-function 
temperature-dependent  conductivity  model  is  unrealistic.  In  fact,  the  local  elec¬ 
trical  conductivity  depends  on  the  degree  of  ionization  a,  as  well  as  temperature, 
so  that  ionization -lags  in  real  nonequilibrium  flows  can  have  considerable  influence 
on  the  variation  of  c  ithin  the  shock  transition. 

The  macroscopic  global  conservation  and  Maxw  i  equations  used  in  the 
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present  analysis  are  developed  in  Sec.  2,  For  an  electrically  conducting  upstream 
state,  these  can  be  integrated  between  upstream  and  downstream  states,  to  give  the 
u.e  al  hydromagnetic  jump  conditions  In  order  to  express  the  dissipation 

fluxes  (i.  e. ,  the  stress  tensor,  heat  flux  vector  and  current  density  vector)  in 
terms  of  lower -order  dependent  variables,  it  is  assumed  first  that  the  electron 
cyclotron  frequency  was  always  much  less  than  the  electron  collision  frequency. 
Secondly,  the  Navier-Stokes  approximation  is  used,  together  with  a  two-temperature 
modification  of  the  Chapman-Enskcg  expressions  for  transport  coefficients.  The 
first  assumption,  which  rules  out  Hall  currents,  is  removable  in  funeral  by^psing 
a  more  general  version  of  Ohm's  law  ,  but  it  is  justified  specificalJy  for  the  flow 
conditions  of  the  calculations  to  be:  presented  later.  It  is  welJ-xnown  that  the  Navier- 
Stokes  appr oximation  is  questionable  in  connection  with  strong  hydrodynamic  shock 
structure  calculations.  Nevertheless,  its  use  in  the  present  study  is  plausible  on  the 
grounds  that  qualitative  misrepresentation  of  the  structure,  of  the  sort  introduced 
by  the  aforementioned  electrical  conductivity  models,  are  unlikely;  moreover, 
Navier-Stokes  equation.,  have  been  used,  with  some  success,  to  study  the  structure 

of  purely  hydromagnetic  shocks  (Marshall19.  Burgers20,  Ludford21,  Germain22’23, 

24  25 

Bleviss  ,  and  Anderson  ). 

In  Sec.  3,  the  ionization  rates  and  transport  properti--  of  partially  ionized 
nonequilibrium  argon  are  developed  in  terms  of  fundamental  collision  cross-sections. 
The  sole  source  of  electrons  is  taken  to  be  collisional  ionization  by  the  reactions: 

kfi  + 

Ar  +  Ar-*  Ar  +  e'  +  Ar,  (la> 


Kfe  + 

e  +  Ar  Ar  +  e  +  e  , 

k  . 
rA 


where  k^,  kf<j  are  the  forward  ionization  rate  coefficients  and  k  k  are  the 

rA  r  c 
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reverse  three-body  recombination  rate  coefficients.  Tlu.  ‘dnetics  of  Eqs.  (la,b) 

were  treated  previously  in  connection  with  flow  in  the  relaxation  zone  of  a  hydro- 
26 

dynamic  shock  .  It  is  assumed  that  the  plasma  remains  quasi-neutral  throughout 
ao  that  elect rogasdynamic  influences  on  shock  structure^  are  negligible  compared 
to  magnetogasdynamic  effects.  Because  the  reaction  rates,  transport  properties  and 
thermodynamics  of  partially  ionized  argon  depend  on  both  electron  and  heavy- particle 
temperatures,  an  appropriate  electron  energy  equation  is  required. 

Sec.  4  deals  with  the  nature  of  the  shock  structure  integral  curves  which  are 
consistent  with  the  gas-ionizing  archetype.  In  this  portion,  it  is  suggested  that  the 
"ZND"  approximation  of  detonation  wave  theory  is  applicable  to  the  present  problem 
and  the  consequences  of  this  representation  are  examined  for  both  "fast"  and  "slew" 
gas-ionizing  shocks.  In  Sec.  5,  selected  numerical  shock  structure  solutions  are 
presented  and  numerical  techniques  are  treated  briefly.  The  conclusions  of  this 
investigation  are  given  in  Sec.  6,  where  the  applicability  and  relevance  of  the 
present  results  are  discussed  and  potentially  profitable  directions  of  future  research 
are  suggested. 
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2.  GLOBAL  HYDROMAGNETIC  EQUATIONS  IN 
OBLIQUE  MAGNETIC  FIELD  SHOCK  GEOMETRIES 


We  shall  be  concerned  here  with  the  distribution  of  flow  variables  within  the 
transition  region  of  the  oblique  gas-ionizing  shock  whose  geometry,  in  shock-frame 
coordinates,  is  shown  In  Fig.  1.  This  shock  may  be  envisaged  as  having  begun  its 
caree'  as  an  ordinary  hydr  odynamic  gas -ionizing  3hock  which  later  "penetrated"  a 
region  of  nonzero  magnetic  .'ield  and  subsequently  attained  a  ste  .dy-state  structure. 
An  (x'.y'.z')  coordinate  system  is  selected  in  which  the  upstream  magnetic  field 
vector  resolves  along  the  x'  and  z'  axes.  For  the  scalar  electrical  conductivity 
assumed  previously,  the  shock-frame  electric  field  is  in  the  y'  direction  and  all 
electromagnetic  components  which  are  initially  zero,  remain  zero  (the  "switch-on" 
shock  is  an  exception  not  explicitly  considered  here). 

As  a  general  rule  "primes"  are  used  to  distinguish  physical  variables,  all 
of  which  are  in  mks  units,  fror-.  the  more  convenient  nondimensionalized  variables 
used  later  in  developing  the  equations;  furthermore,  cartesian  tensor  notation  is 
used  to  express  the  general  form  of  the  governing  equations  more  concisely.  The 
component  directions  in  the  tensor  notation  are  related  to  the  coordinate  system  of 
Fig.  1  as  follows;  x^  =  x',  x'  =  y',  x^  =  z';  (  )T  =  (  )x>  (  >2  =  (  )y,  (  )3  -  (  )z- 
In  the  present  rotation  J.j  is  the  usual  Kroenecker  delta:  1 ,  if  i=j;  if  i/j. 

The  symbol  e.j^  is  the  permucaiion  tensor;  ejsj^O,  if  i= j ,  i=k  or  j=k;  £.^=1.  if  ijk 
are  in  cyclic  order  (123,  32l,  3j  ’)  and  if  ijk  are  unequal  but  not  in  cyclic 

order  (132,  2l3,  32 1) . 

The  thermodynamic  pressure  p'  and  specific  enthalpy  i'  of  partially  ionized 
two  temperature  argon  satisfy  the  equations  of  state^ 


p'  =  C'R(T'  +  aT')  , 


i'=iR(T  +  *Te>  +aR0i'on’ 

where  p'  is  the  mass  density,  T'  is  the  heavy  particle  temperature,  T^ 
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i«  the  electron  temperature,  a  i»  the  degree  of  ionization,  R  =  2.082  x  102 

joule/kg  K  is  the  gas  constant  for  atomic  argon  and  =  183  100  °K 

ion  * 

is  a  characteristic  temperature  for  the  single  ionization  of  argon. 

Using  Eqs.  (2)  and  (3)  to  immediately  eliminate  pressure  and  enthalpy, 
the  global  conservation  and  Maxwell  equations  for  the  steady  flow  of  a  quasi- 
ncutral  plasma  can  be  written  in  divergence  form  as  follows:2^ 


(»'UP  ’  0 


^  1  /  /  5  /  /  ?  .  ,3 

„  I  M  \y  .  _  .  O  /T '  ±  a  r*  '  i  C  ~  '  .  /  /  i: 


7  °Ui  -T'R(T  +  *V  +  f*  a(4n>  +  pV*~ 


E.'B  '  1 

4.  _  Ik.'/,  /  , 

+  +  u.  t..  +  o«  +  q  .  =0 

Uk  u  J  »J  4e,i 

o 


eijk - 7 

ax.' 

J 


*Ek 

Cl->k  bx: 

j 


wo  Ji  * 


Ill— IIWIilllll . 11  « . » 


where  u.  is  the  flow  velocity,  B/  is  the  magnetic  induction,  is  the 

electric  field  intensity,  uq  =  4n  x  10  ^  henry/m  is  the  free  space  magnetic 
permeability  in  mks  units,  is  the  current  flux  density  vector ,  q/  and 
^e  i  are  heavy- particle  and  electron-gas  heat  flux  vectors  and  is 

the  viscous  stress  tensor. 

As  indicated  in  Sec.  1  ,  a  scalar  electrical  conductivity  o  assumed 
in  this  analysis,  in  which  case  the  relevant  form  of  Ohm's  law  is 


Ji  =  a(Ei)  =  °<Ei'  +  «ijkuj'Bk  >  •  <“> 

where  )  =  E.  +  e  .  ^  uj  is  the  electric  field  in  coordinates  moving 

with  the  gas  velocity  u^  through  a  magnetic  field  B^.  Using  the  Navier- 
Stokes  approximation  discussed  in  Sec.  1,  and  recognizing  that  the  partially 
ionized  plasma  is  a  mixture  of  monatornic  heavy  particles  'atoms  and  ions), 
and  an  electron  gas  which  can  in  general  maintain  distinct  temperatures, 
the  heat  fluxes  and  stress  tensor  can  be  written 


qi  = 


qe,i  =  ’  % 


(J& 


n[  +  *1 

l  *x.'  axj 


3  lJ  wj 


wrnmmmsmmm  wmmmmmtmwL  m *—■■«***"-« 


where  x  and  x  are  the  heavy-particle  and  electron- gas  thermal  conductivities 
and  n  is  the  coefficient  of  shear  viscosity  for  the  entire  gas.  Combining 
the  Maxwell  equation  for  induced  magnetic  field,  Eq.  (7),  with  Ohm's  law, 

Eq.  (11),  gives  an  expression  for  the  gas -frame  electric  field  in  term  of 
magnetic  field  derivatives 


Eqs.  (2)  -  (14)  are  applicable  within  the  transition  region  of  Fig.  1.  It 
is  useful  to  re-express  the  governing  equations  in  terms  of  new  "unprimed" 
variables  which  have  been  nondimensionalized  with  respect  to  quantities  in 
front  of  the  shock.  Define; 


-fir  .  T 

Px 


0 


ion  = 


/ 

®ion 


(15) 


'i  * 


Xl 


i 
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where  the  relationship  between  tensor  indices  and  the  components  of  Fig.  1 
has  been  discussed.  Note  also  that  the  subscript  1  in  Eqs.  (15)  -  (18) 
denotes  upstream  conditions  generally,  and  that  i'  the  mean  free  path 
in  the  undisturbed  gas. 

Acoustic  and  Alfv^n  speeds  a^  and  b  are  defined  which  are 
characteristic  of  the  undisturbed  (  =  0  )  state; 
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These,  in  turn,  may  be  used  to  define  the  Mach  and  Alfv4n  numbers  of  the 
shock  Mj^  and  : 


’  Mai 
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In  order  to  assess  the  relative  significance  of  viscosity  verses 
electrical  conductivity  as  dissipative  mechanisms  these  transport  properties 
must  be  incorporated  into  suitable  dimensionless  numbers,  i.e,  ,  fluid  dynamic 
and  magnetic  Reynolds  numbers  R®  and  Rm.  Noting  that  the  characteristic 
length  scale  in  the  present  problem  is  the  upstream  mean  free  path  > 


define; 
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furthermore,  a  Prandtl  number  Pr  is  defined  which  incorporates  the  effects 
of  heavy -particle  thermal  conductivity, 

C  R  M 

Pr  =  4  •  -  •  (22) 

2  K 

It  is  noted  in  passing  that,  from  kinetic  theory,  in  a  pure  monatomic  gas 
v  -  (15  Rfl)  /4,  so  that  when  a  =  0  .  Pr  =  2/3  . 

The  governing  equations,  Eqs.  (4)  -  (14)  can  now  be  written  in  terms 
of  the  dimensionless  quanties  defined  by  Eqs.  (15)  -  (22)  as  follows: 
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Bearing  in  mind  the  relationship  between  the  tensor  notation  indices  and  the 
components  of  the  vector  quantities  in  Fig.  1,  and  substituting  these  quantities 
into  the  above  set  gives  the  ordinary  differential  equations: 


(23) 


d 

d  x 


0  (24) 


Eqs.  (23)  -  (26)  can  be  integrated  immediately  between  conditions  in  the 
undisturbed  gas  and  some  arbitrary  point  x  in  the  shock  interior.  Note 
first  that,  using  the  definitions  of  Eqs.  (15)  and  ( 1 6) ,  the  flow  variables 
must  satisfy  the  following  conditions  asymptotically  upstream 
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Now,  substituting  the  fluxes  of  Eqs.  (28)  ami  (29)  into  Eqs.  (24)  -  (26),  and 
performing  the  aforementioned  integrations  with  the  boundary  conditions  of 
Eq.  (30),  Eqs.  (23)  -  (27)  become: 
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- -  =  Rm  B  (u - L)  +  — ?.i  .  E _ _ L 

dx  L  M.  M,  y  Re  dx 


It  is  instructive  to  examine,  at  this  point,  the  significance  of  the  hydro- 
magnetic  boundary  condition  on  the  electric  field  Ey* .  If  the  boundary  conditions 
of  Eq.  (30)  are  introduced  into  Eqs.  (32)  -  (34),  the  flow  derivatives  quite 
properly  vanish  identically  in  the  undisturbed  gas: 


-  oo  : 


du  du 

_ x  _  _ z_  _  dT 

dx  dx  dx 


In  order  to  insure  that  the  transverse  magnetic  field  vanishes  upstream,  i.e.  , 
@x  - »  -  oo  :  dBz/dx=0,  it  is  required,  from  Eqs.  (30)  and  (35),  that 


Sin  % 


(36) 


@x  -  -co  E.m(B  -E  )  -  0  . 

^  '  zi  yi 

In  a.  pure  hydromagnetic  discontinuity  where  the  gas  is  electrically  conducting 
upstream,  Rm1  j/  0,  so  that  from  Eq.  (36)  Ey^  =  B^  .  On  the  other  hand,  for  the 
gas-ionizing  shocks,  of  interest  here,  q1  =Rm1  =  0,  so  that  E  is  not  uniquely 
defined. 

Eqs ,  (32) -(35)  are  four  differential  equations  in  the  six  unknowns: 
u  .  u  ,  B  ,  T,  T  and  a.  In  order  to  mathematically  close  the  set,  two  additional 
equations  are  required  describing  the  nonequilibrium  behavior  of  a  and  within 
the  shock  transition:  also  the  transport-property-dependent,  dimensionless  numbers 
Re,  Rrr  and  Pr  must  be  expressed  in  terms  of  local  valuer  of  the  flow  variables. 
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3.  IONIZATION  RATES  AND  TRANSPORT 
PROPERTIES  IN  PARTIALLY  IONIZED  ARGON 

Formulation  of  equations  which  deal  specifically  with  distinct  electron, 
atom  and  ion  species  is  facilitated  by  introducing  the  follow  ng  approximations, 
definitions  and  derived  relations,  moat  of  which  follow  directly  from  the 
assumptions  of  Sec,  1; 
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where  n/  ,  n  and  n/  are  the  number  densities  of  Ar  ,  e  and  Ar  species 
respectively,  r_ '  is  the  total  number  density,  nj  is  the  net  electron  number 
density  production  rate  from  all  sources,  j  '  is  the  degree  of  ionization  production 


■31 


-2b 


rate  from  all  sources,  rr>  =  9.107  x  10  kg  and  m.  =  6,628  x  10  kg 

6  A 

are  the  masses  of  an  electron  and  an  argon  itonr.  respective!;  . 

The  one -dimensional  conservation  of  electron  mass  and  energy  equations 
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applicable  to  the  present  oroblem  can  be  written 
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where  k=  1.380  x  10  ^  Joule/ °K  is  Boltzmann’s  constant.  (  nj ).  and 

(n  '  )  are  the  electron  density  production  rates 'resulting  from  atom -catalyzed 
e  o 

reactions  Eq.  (la),  and  electron-catalyzed  reactions ,  Eq.  (lb),  respectively, 
and  v  ’  is  the  collision  frequency  of  the  electron  gas.  The  effects  of  electron 
thermal  conductivity  were  not  included  in  Eq.  (42)  in  anticipation  of  a  future 
development,  however,  a  Joule  heating  term  J  ,  S/  j  was  added  to  the  energy 
equation  of  Ref.  26  to  account  for  dissipation  due  to  induced  currents  flowing 
through  the  gas  within  the  transition  region. 

Making  use  of  Eqs.  (4),  (38)  and  (39).  and  the  fact  that  R=k/m^,  Eqs.  (41) 
and  (4-),  in  terms  of  the  degree  of  ionization  e.  become 
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(43) 
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The  collisional  ionization  source  term.  <„  . 

reactions,  n  ,  and  ^-catalyzed  and  electron. 

A  °e  respectively,  can  be  expree.ed26 
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..  .  y  f  ’  e8re«  of  ionization  a  and  either  h,*  u 

particle  temperature  T'nrtt,  ,  h  heavy- 

P  ature  r  or  the  electron  temperature  -  '  Tu 

v  iure  These. in  turn 
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»e,in  turn,  are  defined  by: 
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thehea15"’11'  ^  K<,q‘T''  *"  "lumbrium  -’constants"  associated  with 
the  heavy- particle  and  electron  th 

n  temperatures,  respectively  Th* 

v  ety.  The  recombination 
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rate  coefficients,  kr  (T  ')  and  k  _  (T  /)  ,  and  equilibrium  constants 
A  ie  “ 

associated. with  Eqs.  (la,b)  behind  strong  normal  shocks,  as  discussed  i.i 
Ref.  26,  are; 


krA<T/>  =  5‘80  X  10~49‘  ‘  exp  (m% 

kre  (Tg)  =  1-29  x  10'44.  ^ii5!3?0  +  2^.  exp^iLi^  (m6/' 


(m~/sec),  (47a) 


sec)  ,  (47b) 


Kep  (T ')  =  2.90  x  1022.  T '  ^3  *  exp  (-  (l/m3)  , 


Keq  (Tj)  =  2.90  x  1022  •  Tj’4  • 


(l/m3). 


It  is  convienient  to  define  nondim ensionalized  "unprimed"  variables 
"A  ’  ae  ’  v  »  corresponding  to  the  production  rates  and  collision  frequencies 
appearing  in  Eqs.  (43)  and  (44): 
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The  production  rates  of  Eq.  (49)  are  completely  specified  in  terms  of  local 
values  rf  u  ,  T  and  T  *  by  Eqs.  (45a)  though  (48b).  Relations  are  now  sought 
which  express  and  also  the  dimensionless  numbers  Re,  Rm  and  Pr  in 
terms  of  a  .  T'and  T  '  . 

In  principle,  all  collision-dependent  transport  properties  needed  in 
this  analysis  are  obtainable  from  a  knowledge  of  the  elastic  collision  cross 
section  for  the  various  encounters  occuring  in  a  partially  ionized  gas.  These 
will  be  briefly  summarized  for  argon. 

The  Coulomb  cross-sections  for  collisions  between  charged  particles 
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where  and  are  the  cross-sections  for  Ar+-Ar+,  e~  -  Ar  + 

and  e  -  e  collisions  respectively,  e=  1,602  x  10-1^  coulomb  is  the 
charge  of  an  electron  and  Pq=  8.854  x  10‘12  farad/m  is  the  dielectric 
permittivity  of  free  space  in  mks  units.  Making  the  required  numerical 
substitutions  in  the  above  yields 


1.95  x  10 


-10 


1.53  x  10 


14  T 


, '  » 


(m*),  (51a) 


M 


■^el  =  3ee 


1.95  x  10 


-10 


'  a 


.  i  n  [l .  53  x  I014.i-] 


(m3)  ; 


(51b) 


the  remaining  cross-sections  can  be  expressed,  after  Jaffrin^; 
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where  and  are  the  elastic  cross  -section  for  Ar  -  Ar  , 

Ar*  -  Ar  and  e  -  Ar  collisions,  respectively. 

The  electron  elastic  collision  frequency  \>J  and  the  electrical 
conductivity  of  the  partially  ionized  gas  o  can  be  written  directly  in  terms 
of  these  cross  sections 
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moreover,  the  thermal  conductivities  of  the  atom,  ion  and  electron  species 
30 
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The  viscosity  coefficients  are  related  to  the  thermal  conductivities  by 
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The  upstream  mean  free  path  X  ^  which  is  used  here  as  a  reference  length  is 
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where  0^1  is  Eq.  (5lc)  evaluated  at  T  '  =  T^  .  Note  also  that 


*a  +  "i  ,  n  =  nA  +  ^  ; 


(59) 
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where  the  electron  viscosity  has  been  dropped  from  ry  since  ry  /  ry  is  of 
L  * 

the  order  (m  /m.)  “  1  ;  cf.  Eqs.  (56)  and  (57). 

6  A 

Combining  Eqs,  (20),  (39),  (54),  (55),  (56),  (58)  and  (59)  and  i-  troducing 
the  definition  of  the  upstream  Mach  number  =  u  ;  /  (5  RTt  )  ~  gives 

an  expression  for  the  fluid  dynamic  Reynold's  number  in  terms  of  a  and  T: 


R  =  -Hi)  M,  T-‘  .  J_  _ Ll« _  +  _ 2 - 

V5"/  QaX  +“(Qlk-QAi>  Qll+  ~-'QlI-QIA 


Using  the  fact  that  =  (T'/t')  ‘  =  T"  **  ,  from  Eq.  (5lc), 

Eq.  (60)  can  be  simplified  cons.  *ably  if  the  argon  remains  un-ionized; 


a  =  0  :  R  =  (iii \ 

e  v5  7 


MjT:3/4 


Combining  Eqs.  (2i),  (39),  (53)  and  (58),  and  using  the  Mach  number,  as 
before,  gives  an  expression  for  magnetic  Reynolds  number  in  terms  of  a  and 


Rm  =  p  M  x  a  T 
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where  p  is  a  nondimensionalized  reference  density  defined  by 
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It  may  be  noted  here  that  while  the  fluid  dynamic  Reynold's  number 
bared  on  mean  free  path  Re  is  density  -  independent,  this  is  not  the  case  for 
the  Magnetic  Reynolds  number  Rm.  From  Eq.  (6l)  it  follows  that  p  and 
therefore  Rm  increase  with  decreasing  upstream  density  , 

Combining  Eqs.  (39),  (50)  (52)  and  (58)  with  the  definition  of  Mach 
number  gives  the  following  expression  for  nondimensionalized  electron  collision 
frequency  in  terms  of  a  and  T  : 


QeA+  °<QeI  •  QeA> 


.  (62) 


From  Eqs,  (22),-  (57)  and  (59)  the  Prandtl  number  of  the  Mixture  is  simply 


pr  =  — - 

3 


The  conservation  of  mass  and  energy  equations  for  the  electron  gas 
can  be  put  in  a  dimensionless  form  consisted  with  that  of  the  global 
conservation  equations  of  Eqs.  (32)  -  (35),  Note  first  from  Eqs.  (17),  (21) 
and  (27).  That  the  term  J  3/a  can  be  written 


_  _  /  p  / 

_  >  ?  RmI3  u 
J  Xi  xi 


( u  B  -  u  B  -  E  ) 
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Substituting  Eq.  (64),  together  with  the  dimensionless  production  rates  and 
collision  frequencies  of  Eqs.  (49)  and  (50)  and  the  unprimed  variables  defined 
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earlier  by  Eqs.  (15)  -  (19),  into  Eqs.  (43)  and  (44)  yields32 

d  a  _  aA  *  “e 
d  x  u 

x 


(65) 


It  is  significant  that  Eq.  (66)  is  not  singular  when  q  =  0  since  from  Eqs.  (45b), 
(49),  (61)  and  (62)  it  follows  that  ve,(ae/o)  and  (Rm/a)  are  all  bounded 
as  r- - *  0.  Certain  formulations  of  the  electron  energy  equation  which 

Live  appeared  in  the  literature  have  not  had  this  useful  and  physically  reasonable 
property. 
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4.  INTEGRAL  CURVES  IN  THE  "2ND"  lPPROXIMATION 


Equations  (32)  -  (35),  (65)  and  (66),  together  with  the  auxiliary 

algebraic  expressions  for  Re,  Rm  ,  Pr ,  a,  ,  and  v  developed  in 

the  preceeding  section,  form  a  mathematically  closed  set  of  six  ordinary 

differential  equations  in  the  six  primary  variables  ux  >  uz  >  T  ,  Tg  and 

e  .  The  formal  solution  of  these  equations  as  an  initial-value  problem  starting 

from  the  upstream  boundary  conditions  of  Eq.  (30)  is  not  possible  however 

owii.^  to  the  mathematical  nature  of  the  system.  Briefly,  this  can  be  explained 

as  follows.  The  leading  edge  of  the  gas-ionizing  front  must  begin  as  an 

ordinary  hydrodynamic  shock,  but  a  predominant  characteristic  of  the  Navier- 

Stokes  hydrodynamic  shock  structure  in  a  monatomic  (Pr  =  2/3  )  gas  is 

that  the  integral  curve  solution  in  (ux  ,  T)  phase  space  has  a  singularity  of 

the  node  type  at  the  upstream  state;  consequently,  the  downstream  state  is 

33 

"unattainable"  from  the  upstream  state  by  numerical  integration  .  Since 
the  gas-ionizing  shock  begins  its  upstream  structural  development  as  a  pure 
gasdynamic  shock,  the  latter  conclusion  applies  to  the  present  case  as  well. 

Fortunately.it  is  appropriate  to  employ  a  useful  approximation  here  which 
has  been  developed  in  the  theory  of  detonation  waves.  Commonly  known  as 
the  Zeldovich-von  Neumann-Doring  (ZND)  approximation,  in  the  present  context 
this  amounts  to  recognizing  that  ionizing  reactions  of  Eqs.  (la,b)  are  sufficiently 
"slow"  such  that  gas -ionizing  shock  structure  can  be  computed  in  two  distinct 
regions:  (1)  a  perfect-gas  viscous  shock  wave  standing  in  fro”?  of  (2)  a  much 
longer  ionization  relaxation  zone  where  finite- rate  chemistry  and  hydromagnetic 
interactions  are  significant. 

22  24  18 

It  has  been  recognized,  e.g.,  by  Germain  ,  Bleviss  and  Leonard  , 

that,  for  hydromagnetic  shocks,  when  the  magnetic  Reynold's  number  is  small 
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compared  to  the  viscous  Reynolds  number,  a  viscous  shock  is  imbedded  in  a  much 

v/ide r  region  of  hydrom.  netic  interaction.  For  the  (collisionally  ionizing)  gas- 

ionizing  hydromagnetic  vaves  treated  in  the  present  paper,  this  must  be  the  case 

since  realistic  ionization  rate-processes  yield  values  of  Rin/Re  «  1  within  the 

initiating  perfect-gas  shock  regardless  of  the  ultimate  electrical  conductivity 
34 

level  .  It  should  also  be  clearly  understood  t’-at  unlike  certain  imbedded  viscous 
shocks  which  can  occur  in  pure  hydromagnetic  wave  fronts,  the  imbedded  shock 
here  must  stand  upstream  of  the  hydromagnetic  interaction  since  it  creates  the 
necessary  electrically  conducting  environment. 

The  equations  governing  flow  in  the  two  regions  c  n  be  obtained  formally 
from  Eqs.  (32)-(35),  (65)  and  (66)  by  applying  the  appropriate  limiting  conditions. 

In  the  perfect-gas-shock  region  we  have  the  limit:  a  -  0,  Rm  -  0,  so  that  Eqs.  (65) 
and  (66)  for  a  and  Tg  are  not  relevant  a.id  Eqs.  (30),  (32) -(35)  become 


du 


dx 


=  0  ,  u  =  u  =  0  ; 
*  z  zi 


dB 

z 

dx 


=  0  ,  B  =  B 
*  z  zi 


(66) 


du 


x 


dx 


(67) 


dT 

dx 


Re 


|T-1  + -  (u. 


(68) 
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where  we  have  used  ths  fact  that  Pr  =  2/3  in  the  above.  Clean  /  ihe 

transverse  magnetic  field  and  transverse  velocity  are  constant  across  the 

perfect-gas  shock.  Eqs .  (67)  -  (68)  describe  the  ^vier-Stokes  shock 

structure  of  a  perfect  monatonuc  gas.  Their  solution  has  been  treated  eloe- 
35-37 

where  ,  and  is  discussed  here  in  the  Appendix.  The  upstream  and  downstream 

states  implied  by  these  equations  can  be  found  by  setting  du^  /  dx  =  dT/dx  =  0 
in  Eqs.  (67)  and  (68), 


3  T  m 

u  -  1  +  - —  ( - l)  =  0  .  T  -  1  +  — i.  (u*  -  1)  =  0, 

5M3  u  3  x 

l  x 


eliminating  the  temperature  T  between  these  to  get  the  quadratic 


(*'=r)  ■■’('*  =r) ' 


and  solving  for  the  velocities  and  corresponding  temperatures  associated  with 
the  two  roots.  This  yields  the  upstream  and  downstream  states  of  the  imbedded 
parfect-tas  shock: 


-  oo  ;  ux  =  T  =  1  , 


I  u  =  -  (l  +  .  T  =  -1-  (M3  -  -L)  +  1  . 

x  4  vfa  16  x  ct*a'  8 


5M"  8 


In  the  relaxation  zone  regions  the  governing  equations  are  found  by 
applying  the  limit  Re — -*  oo  to  Eqs.  (32)  -  (35)  which  yields  the  set 
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T  +  aT 


ux  -  1  +  ~ 
x  5M8 
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-  1  + 


2M 


(b!  -Bj)  =  0, 


3  Z 


(70) 


a  s 
z 


B  -  B  , 

7  Zl 


(71) 
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Ai 


T  +  ftT.  +  ■=  a 


e  5  ion 


M' 

1  +  —2 -  (U*  -  1) 

3  x 


dB, 

dx 


-  Rm 


+  2  .  K_ 

3  MAi 


B  (u  -  JL 
2  l  X  w  • 

mAi; 


(Bz  -B,lMEy  *  —  B--“  Is0'  <72> 


2M 


Al 


B 


_iL  -  e 

»A;  yl 


=  Rm.g(ux,  Bz)  (73) 


Eliminating  the  quantity  (T  +  aTg)  between  Eqs.  (70)  and  (72)  gives  a 
quadratic  equation  in  u^  corresponding  to  Re  — >  co  , 


f(ux  •  Bz  ’  ^  4ux  -(5  +  +  £1  lux+  1  +  .-4-  -  e  I  =  0*  <74> 


M' 


M 


where 


5 

2 


Ba  -  B  8 
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M 
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—  .  - —  (B  -  B  )  E  + 
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B  -  B 

Z  21 


In  the  ionization  relaxation  zone  the  flow  will  proceed  along  the  path  f  =  0; 
it  is  useful  to  introduce  an  analogous  path,  g  =  0,  corresponding  to  the  limit 


Rm-*0 


g  (  u  B  )  =  B  [  u  - 

11  2  V  *a\ 


— -  -  E  =  0. 

w  3  yi 

mai 


Equation  (74)  has  two  roots,  given  by 


u  =  J-|5  +  _L_  +  Gl  + 
8  M,s 


+  16  e 


while  solving  for  u  from  Eq.  (75)  gives 


Prior  to  discussing  numerical  shock  structure  solutions  in  physical 

space,  it  is  instructive  to  examine  the  path  of  the  ZND  solutions  in  (u  B  1 

'  x  ’  z' 

phase  space.  Since  the  integral  curves  are  more  mean’ngful  if  a  distinction  is 
made  between  "fast"  and  "slow"  hydromagnetic  waves,  these  classifications  will 
be  briefly  reviewed  in  the  context  of  the  present  work. 

Analysis  of  the  linearized  hydromagnetic  equations  yields  the  so-called 
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fast  and  slow  disturbance  speeds  Cj  '  and  ca  which,  together  with  the 
acoustic  and  Alfven  speeds  a^  and  bXj  of  Eq.  (19a),  are  properties  of  the 
undisturbed  flow.  These  speeds  are  conveniently  written  in  terms  of  the 


quantity 


■  i  /  \  g  ■ 

b'  _  (b  '*  +  b  ,9)>5  =  b  '  1  +  (-^-1  =  b  ' ,  (1  +  B  s  )** 

1  =  '  xi  zi  '  xi  I  g  il  xi  '  zr 

\  xi/  - 


as  follows:" 


'  1  1  '  a 


/a,,  /a, a  .  /a.  >a 


f  -  i  (aT+bi-)  +>/(ai-+lhT- 


4a,  b 


'  1  /a,  ,  /a,  j,  /a  ,  ,  >a.a  .  'a.  >  ■ 

j=  —  (at  +bj  )  -  V (ai  *  )  *  ^ai  1>X1 


/a  ,  v'a.a  A<  a  ,  ’  a 


Pure  hydromagnetic  shocks  are  generally  classified  as  either  fast  or  slow 

39 

depending  on  whether  they  satisfy  the  inequalities 


xi  >  l  ;  fast  shock, 


<  1  ;  slow  shock. 


Using  ‘he  definitions  5  u^  /  a^  ,  -  u^  AX  l  and  Eq,  (78), 

the  following  useful  formulas  are  obtained; 
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For  the  purposes  of  this  section,  attention  is  restricted  to  (hydromagnetically) 
oblique  shocks  in  the  infinite  Mach  number  limit,  since  hypersonic  (Mi»  1)  Mach 
numbers  yield  ionization  levels  required  for  hydromagnetic  interaction  and  this 
particular  limit  does  not  change  any  important  features  of  the  integral  curves.  For 
Mx  — •  oo  then,  Eqs.  (81a, b)  give  u^/c^  M^I  +  bM"1/3  and  Csi/Uxi  =  °*  For 
the  special  case  when  the  upstream  magnetic  field  is  inclined  at  45°  to  the  shock 
front  (B  =  B  =  1),  the  criteria  of  Eqs.  (80a,  b)  become 

Xl  T»  1 


M^  >  ^"7  :  fast  ohock, 

0  <  M.  <  1  :  slow  shock. 

—  Ai 


(82) 

% 

§ 

t 


It  might  be  observed  here  that  no  pure  hydromagnetic  shocks  can  exist  between 

the  weakest  (acoustic)  fast  wave  at  and  the  slowest  (switch- off) 

slow  wave  at  M.  =  1 . 

A  i 

Returning  to  the  discussion  of  integral  curves  in  the  ZND  approximation, 
the  45°  upstream  magnetic  field  and  M^  — »  oo  assumptions  (which  were  introduced 
to  make  the  problem  specific)  should  be  borne  in  minu  as  they  apply  to  the  balance 
of  this  section.  From  Eqs.  (69a,  b)  one  can  expect 
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an  initial  jump  in  streamwise  velocity  from  u^  =  i  to  =  1/4  whiie 

B  =  1  remains  constant,  corresponding  to  the  perfect-gas -shock  transition, 

z 

$  $ 

Tlris  imbedded  shock  transition  is  denoted:  1  -  1  ,  where  1  is  the  downstream 
state  of  the  perfect-gas  shock.  Subsequently,  the  flow  progresses  in  the 
relaxation  zone  along  the  path  f  (u^  ,  B  ,  a)  =  0  ,  until  the  downstream 
state  of  the  gas -ionizing  shock  is  attained.  The  latter  step  of  the  overall  shock 

j ;« 

transition  is  denoted  ;  i  -  2.  As  a  conseqence  of  all  flow  derivati  es  vanishing 

downstream,  e.g.  du  /  dx  =  dB  /dx  =  0  @x - )  oo,  the  downstream  state 

in  the  (u^,  B^)  plane  is  indicated  by  the  intersection  of  the  curves  f  (u^  ,  B^)  =  0 

and  a  (u  ,  B  )  =  0.  Note  that  in  the  infinite  Mach  number  case 
6  '  x  z 

6  o0. 

lim  M  - 1  oo  - —  =  0  , 

5  M“ 


so  that,  from  Eq.  (74),  f=f(u  ,B  )  does  not  depend  on  a. 

X  z 


(a)  Fast  Shocks 

Consider  now  the  possible  trajectories  of  the  fast  shock  =  J  l  O'  (>^2) 

shown  in  Fig.  2.  Recall  from  Eq.  (36)  that  because  **  =  Rm  =  0  ,  an 

i  i 

ndeterminacv  exists  in  the  value  of  E  for  gas-ionizing  shocks;  it  is  there- 

y  i 

ore  appropriate  at  this  point  to  treat  the  shock-frame  electric  field  as  a  free 

arameter.  This  was  done  in  Fig.  2  which  shows  the  curves  f  =  0  and  g  =  0 

omputed  from  Eqs .  (76)  and  (77)  for  (a)  E^  =  1.0,  (b)  E^  =  0.625  and 

c)  E  =  0.25.  The  first  trajectory,  nam  that  with  E  =  B  =  1.0, 
yi  yizi 

orresponds  to  the  upstream  boundary  condition  on  thn  electric  field  in  a  pure 
ydromagnetic  shock,  cf.  Eq.  (36).  As  indicated  previously,  the  transition,  if 
;  occurs,  must  take  place  by  the  path  1-1  -  2  in  the  ZND  model.  It  can  be 
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shown  that  this  path  is  impossible  for  fast  shocks  from  the  following  argument; 

Since  g  <’  0  below  the  curve  g  =  0,  it  follows  from  Eq.  (73)  that  dB  /dx  <r 

z 

>Js 

0  at  point  1  (Rm  >  0,  of  course)  ;  but  >  1  ,  from  the  intersection  point 

of  g  =  0  and  f  =  0  ;  therefore  ,  the  downstream  state  is  inaccesible  by  the 

*!< 

path  1  -  2  along  the  f  =  0  curve  since  the  magnetic  induction  equation 

predicts  a  decrease  rather  than  the  required  increase  in  transverse  magnetic 
field.  The  same  reasoning  applies  for  all  valves  of  >  1/4  ,  cf.  Fig. 

2(b). 

An  analogous,  but  oppositely  directed  situation  occurs  when  E^  <r  1/4, 
since  B^  <  1  and  the  1-1-2  transition  becomes  impossible  because 

jjc 

g  >  0  along  1  -  2  and  Eq.  (73)  predicts  an  increase  of  B  instead  of  the 

z 

required  decrease.  In  fact,  the  only  permitted  1-1-2  transition  in  a 
fast  gas-ionizing  shock  is  the  degenerate  case  of  E^  =  0.25  in  Fig.  2(c), 
which  is  nothing  more  than  a  pure  gas  shock  with  no  change  in  magnetic  field; 

1). 

If  the  gas  were  electrically  cor-,"rt»-..?  -  stream,  but  with  Rm/Re  <  <  1  , 

the  pure  hydromagnetic  transition  i  -  2  ,  along  the  upper  branch  of  the 

f  =  0  curve  in  Fig.  2(a),  would  be  indicated.  In  this  regime  g  >  0, 

dB  /dx  >  0,  and  there  are  no  contradictions  of  the  type  encountered  in  the  ZND 
z 

gas-ionizing  integral  curves.  As  indicated  previously,  this  branch  must  be  ruled 
out  here  since  it  violates  the  gas-ionizing  archetype  of  Sec.  1. 

Although  switch-on  (B^=0,  B^/0)  and  transverse  (Bx=Bx^=0)  gas -ionizing 
shocks  have  not  been  dealt  with  specifically,  it  can  be  shown  that  ZND  structures 
are  impossible  in  these  shocks,  for  the  same  general  reasons  that  were  given  for 
the  fast  oblique  shocks  discussed  in  this  section. 


(u  =  1/4  ,  B  =  B 

'  X3  '  Zs  Zl 


m 
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(b)  Slow  Shocks 


From  Fig.  3  it  is  evident  that  the  f  =  0  and  g  =  0  curves  for  slow 
shocks  take  substantially  different  forms  from  those  of  the  aforementioned 
fast  shocks.  These  curves  are  plotted  for  the  slow  shock  i//T  (<d 

at  five  difierent  values  of  the  shock-frame  electric  field:  (a)  E  =0.25, 

y  i 

(b)  E  =  0.625.  (c)  E  =  1.0,  (d)  E  =  1.50  and  (e)  E  =  2.0. 
yi  yi  yi  yi 

In  this  case,  Hg.  3(c)  is  the  plot  associated  with  the  pure  hydromagnetic 
boundary  condition  on  the  electric  field.  All  the  electric  fields  shown  have 
in  common  the  property  that  ZND  structures  are  possible.  After  the  1  -  1 
perfect-gas-shock  transition  the  flow  is  in  a  region  where  g  <  0  ;  but  Bzg  <r  1 

in  this  case  so  the  derivative  dBz/dx  <^0  [from  Eq.  (73)]  is  in  the  proper 

$ 

direction,  thereby  permitting  the  1-1-2  path  to  the  downstream  state.  There 

are. two  limiting  cases  of  interest:  (1)  E  =  l/4  in  Fig.  3(a)  is 

the  degenerate  case  corresponding  to  the  pure  gasdynamic  shock  with  constant 

magnetic  field  and  (2)  E  =  2.0  in  Fig.  3  (e),  which  yields  a  new  kind  of 

switch- off  shock  which  can  only  occur  in  gas -ionizing  fronts.  It  follows  from 

Eq.  (73)  that  when  E  =  B^  /  =  2.0,  g  =  0  along  the  straight  lines 

B^  =  0  and  u^  =  1  /  =  2  .  The  latter  part  of  the  g  =  0  curve  is  not  visible 

in  this  plot  because  the  ordinate  is  cut  off  at  u  =  1.2  .  Since  the  intersection 

of  the  f  =  0  curve  with  B  =0  corresponds  to  the  downstream  stata,  the  lab- 

z 

fiame  electric  field  E  =  B  /  M.  3  results  in  a  complete  switch-off  of  the 

yi  zi  '  Ai  r 

transverse  magnetic  field.  The  sv/itch-off  shock  is  only  possible  in  ordinary 
hydromagnetics  when  =  1,  but  it  is  clearly  obtainable  in  slow  gas-ionizing 

shock  waves  propagating  at  other  Alfv^n  numbers. 
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5.  NONEQUILIBRIUM  NUMERICAL  SOLUTIONS 


It  should  be  clear  by  this  point  that  the  mathematical  nature  of  the  problem 
in  the  present  ZND  approximation  is  fundamentally  different  from  that  implied 
originally  by  Eqs.  (32) -(35),  (65),  (66)  and  the  associated  initial  conditions  of 
Eq.  (30).  Rather  than  attempt  the  solution  of  six  differential  equations  in  the  major 
flow  variaoles  u^,  u^,  Bz ,  T,  Tg  and  a,  it  la  proposed  instead  to  solve  simpler  sets 
of  equations  in  two  different  regions  and  to  match  their  solutions  at  a  suitable  point. 
Specifically,  in  the  perfect-gas -shock  region,  only  two  differential  equations  need 
be  integrated  [Eqs.  (67)  and  (68)],  while  In  the  ionization  relaxation  zone  there  are 
three  [Eqs.  (65),  (66)  and  (73)]  ,  In  the  latter  case,  the  velocity  components  u^,  uz 
and  the  heavy-particle  temperature  T  are  evaluated  locally  from  algebraic  relations 
[Eqs.  (71),  (72)  and  (76)].  In  addition,  the  local  values  of  ct^,  'nd  v  needed 

to  numerically  integrate  the  relaxation  zone  differential  equations  are  available  from 
relations  introduced  and  developed  in  Sec,  3  [Eqs.  (45a, b),  (51b, c,d),  (6la,b)  and 
(62)]  . 

In  order  to  solve  for  the  distribution  of  flow  variables  within  the  gas -ionizing 
front,  two  different  IBM  7090  computer  programs  were  created;  one  to  solve  the 
perfect-gas  Navier-Stokes  shock  structure  problem  (see  Appendix)  and  the  other  to 
solve  the  hydromagnetic  ionization  relaxation  zone  problem.  Since  the  perfect-gas  - 
shock  structure  extends  from  x=-ootoxs+oo  and  the  ionization  relaxation  extends 
from  some  finite  value  of  x  (say  x=0)  to  x*  +  oo,  the  two  regimes  overlap  in  physical 
space;  consequently,  it  was  necessary  to  cut  the  perfect  gas  shock  solution  off  at 
some  arbitrary  point,  as  explained  in  the  Appendix  and  tack  it  on  again  to  the 
beginning  of  the  relaxation  zone  in  order  to  construct  a  single-valued  solution 


over  the  entire  range  of  x,  from  ~ao  to  +  oo.  The  somewhat  legislated  nature  of  this 

matching  procedure  is  characteristic  of  so-called  singular  perturbation  problems  — 

these  are  invariably  generated  by  physical  processes  which  involve  disparate  length 

scales  —  and  is  a  familiar  feature  of  boundary- layer  solutions  and  detonation  wave 

structure  solutions  in  the  ZND  approximation  .  It  is  possible,  in  principle,  to 

obtain  a  more  rigorous  formulation  of  the  connecting  region  between  the  two  solutions 
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by  the  method  of  matched  asymptotic  expansions 

For  all  calculate  i  discussed  in  this  section  it  was  assumed  that  the  fronts 
propagate  into  "cold"  un-ionized  argon  which  is  at  a  pressure  of  pi  =  1.0  mm  Hg  = 

1.33  x  102  newton/ma  and  a  temperature  of  Ti  =  300°K  and  has  a  corresponding 

_5 

upstre  nn  mean  free  path  of  =  5.  38  x  10  mj  furthermore,  all  of  these  fronts 
were  considered  to  be  traveling  at  the  same  gas  dynamic  Mach  number  Mi  =20, 
cor  responding  to  an  upstream  flow  velocity  of  6.45  x  103m/sec,  when  viewed  from 

a  shock-fixed  reference  frame.  It  will  become  evident  that  a  number  of  phenomena 
of  interest  develop  at  these  flow  conditions. 

Figure  4  shows  the  results  of  a  combined  perfect-gas-shock  and  relaxation 
zone  calculation  for  the  limiting  case  of  an  ordinary  hydrodynamic  gas-ionizing 
front,  i.e.  ,  with  no  imposed  electric  or  magnetic  fields:  E  =  B  =  B  =  0. 

The  solutions  have  been  jointed  at  x  =  0  and  the  scale  has  been  stretched  by  a 
factor  of  ten  for  x  <  0  compared  to  x  >  0  scale  in  order  to  show  the  relatively  narrow 
Navier-Stokes  shock  structure.  Clearly,  the  ionization  relaxation  takes  place  over 
several  hundred  upstream  mean  free  paths  compared  to  the  few  mean  free 
paths  required  by  the  perfect-gas  shock,  thus  providing  an  L  posteriori  verifica¬ 
tion  of  the  assumptions  leading  to  the  ZND  approximation.  The  relaxation 

26 

zone  behavior  of  a  purely  hydrodynamic  front  has  been  discussed  elsewhere 
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It  should  be  mentioned  that  the  initia1  electron  temperature  used 

here  in  the  numerical  integration  of  Eq,  (66)  in  all  cases  was  taken  as  (0)  = 

:'c 

T  (0)  =  T  .  Although  the  value  of  electron  temperature  immediately  behind 
the  perfect-gas  shock  is  not  well  defined  (since  n  =  0  there),  it  was  shown  in 
Ref.  26  that  relaxation  zone  calculations  are  almost  entirely  insensitive  to 
arbitrarily  selected  initial  values  of  the  electron  temperature. 

For  reasons  explained  in  Sec.  4,  it  is  not  possible  to  compute  ZND 
structures  for  fast  gas -ionizing  fronts,  so  attention  has  been  turned  toward 
the  family  of  slow  oblique  shocks  discussed  previously,  whose  upstream 
magnetic  field  is  inclined  at  45°  to  the  front  and  whose  Alften  number  is  = 

1  /  J  2.' .  In  dimensional  terms,  the  corresponding  streamwise  magnetic 
field  upstream  is  B  =  0.473  Wb/m^,  a  reasonably  attainable  value  in 
laboratory  experiments.  Structure  calculations  were  carried  out  with  various 
values  of  the  electric  field.  These  results  are  displayed  in  Figs.  5  and  6, 
respectively,  for  the  pure  hydromagnetic  boundary  condition  on  the  electric 
field  =  1.0,  and  for  the  gas-ionizing  switch-off  shock  E^  =  2.0,  whose 

■omewhat  unique  existence  was  discussed  earlier  in  Sec.  4.  As  the 
nonequilibrium  ionization  progresses  and  the  gas  becor.  es  electrically  conducting 
the  transverse  magnetic  field  is  decreased.  In  the  case  of  Fig.  6  it  is  completely 
switched  off.  The  energy  associated  with  the  magnetic  field  is  consequently 
transferred  into  other  modes,  i.e.  ,  thermal  and  ionization  energy.  The  cc- verted 
magnetic  energy  can  be  viewed  as  an  effective  exothermicity  within  the  front. 

The  scale  stretching  for  x  <  0  discussed  previously  was  also  applied  to  the 
plots  in  Fig.  5  and  6  so  that,  even  though  they  appear  to  be  approximately 
the  same  width,  the  perfect-gas -shock  is  still  an  order  of  magnitude  narrower 

than  the  relaxation  zone  in  the  extreme  case  of  E  =  2,  0  Fig.  6  ). 

Yi  8 
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As  indicated  in  Sec.  1,  it  has  been  held  by  certain  writers  that  the  electric 


field  Ey  associated  with  a  gas-ionizing  hydromagnetic  front  will  be  determined  by 
the  structure.  (  r.  the  other  hand,  it  was  shown  in  Sec.  4  that,  although  structural 
considerations  may  well  rule  out  the  steady-state  existence  of  certain  (fasti  gas- 
ionizirg  shocks,  they  do  not  appear  to  furnish  a  criteria  as  to  which  of  the  possible 
electric  fields  will  actually  be  observed  in  the  (slow)  shocks  whose  existence,  in 
^  ZND  approximation,  is  possible.  Fig,  7  illustrates  computed  distribution  of 
downstream  /alues  ,  u^  and  a3  corresponding  to  the  upstream  conditions 


XB 


discussed  previously  for  various  values  of  ‘he  shock-frame  electric  field  .  Also 

shown  is  a  scale  indicating  the  corresponding  nondimensionaiized  lab-framt  >  Lectric 

field  (which  happens  to  equal  the  upstream  gas-frame  electric  field  Eyi=  E  -B^, 

since  the  undisturbed  gas  is  obviously  motionl  ss  with  respect  to  the  laboratory). 

Ev’dently  Bzs  decreases  almost  linearly  with  increasing  electric  field  until  it  is 

finally  switch'd- off  at  E  =2.0.  Tt  is  interesting  that  B  =1.0  occurs  at  E  ~  0. 1, 

Yi  zs  yi 

rather  than  E^=0.25  as  one  might  expect  from  Fig.  3(a).  This  is  due  to  the  finite 

Mach  number  used  in  the  present  calculations,  so  that  6  a®.  /5M  *  was  not  zero 

'e  assumed  in  the  Fig.  3  plots.  As  decreases  in  Fig.  7,  aa  and  Tg  increase  as 

energy  is  redistributed.  Ultimately  when  the  equilibrium  gas  becomes  fully  ionized, 

at  about  E  =1.2,  the  kinetic  energy  u  a /2  and  hence  u  increases  as  well, 
yi  *  x'  x 

The  extent  of  the  relaxation  zone  can  be  estimated  if  we  define  a  suitable 
characteristic  length  /  (strictly  speaking,  of  course,  equilibrium  is  not  attained 
until  x  -  oo).  Consistent  with  Ref.  26,  let: 


l  ?  [xl  L 

r  Uo=(?)  * 


(83) 


The  upper  graph  in  Fig.  7  shows  the  variation  of  this  nondimens ioualized  relaxation 
length  ,  where  i  '  is  the  physical  relaxation  length.  Since  the  perfect-gas 

shock  is  the  same  thickness,  <-s  ts/\  =  7,  in  all  cases  (because  it  depends  only  on 
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6.  CONCLUDING  REMARKS 


The  present  study  has  dealt  theoretically  with  the  existence  and  structure  of 
gas -ionizing  hydromagnetic  shock  waves ,  as  defined  by  the  archetype  of  Sec.  1. 

In  the  case  singled  out  for  attention,  the  gag'  within  the  front  was  collisionally 
ionizing  argon  (by  atom -atom  and  electron-atom  impacts)  in  a  nonequilibrium  two- 
temperature  state.  In  view  of  the  relatively  low  temperature  expected,  photo- 
ionization  was  ruled  out  on  an  ad  hoc  basis;  moreover,  the  ionization  lags  associated 
with  finite- rate  chemistry  indicated  that  the  ZND  approximation  could  be  employed. 
This  model,  in  turn,  led  to  a  number  of  surprising  results;  (1)  No  steady-state 
structure  could  be  constructed  for  fast  gas-ionizing  waves;  (2)  for  the  slow  waves, 
where  numerical  solutions  were  obtained,  the  ZND  approximation  was  verified 
a  posteriori  fur  the  shock  conditions  studied  here;  and  (3)  the  Rankine -Hu,;oniot 
indeterminacy  o»  the  electric  field,  which  is  intrinsic  to  the  concept  of  gas -ionizing 
shocks,  was  not  removed  by  considerations  of  structure. 

As  to  th>  applicability  and  relevance  of  these  results,  it  should  be  first 
recalled  that  "high  temperature"  gas -ionizing  shocks  were  ruled  out  at  the  outset 
on  the  physical  grounds  that  they  create  radiation-induced  electron  precursors  and 
hence  make  the  upstream  state  electrically  conducting;  but  any  "gas -ionizing"  shock 
which  propagates  in  the  real  world  must  move  into  a  region  where  there  is  some 
elec  *ical  conductivity,  however  small.  Consequently,  the  gas-ionizing  archetype 
actually  presupposes  some  low  threshold  below  which  the  gas  acts  as  though  it  were 
I'.on-electricaliy-conducting.  It  is  not  entirely  clear  that  such  a  threshold  exists 
physically.  Furthermore,  even  in  tig  supposed  low  temperature  case  considered 
here,  the  post  shock  temperatures  become  sufficiently  high  (since  the  slow'  shock 
is  effectively  exothermic)  to  indicate  a  relatively  high  level  of  precursor  ionization, 
particularly  for  large  electric  fields.  For  this  reason,  in  lieu  of  specific  experi¬ 
ment  evidence  to  the  contrary,  it  seems  quite  possible  that  soJutims  with  arbitiarily 
selected  electric  fidefe  aic  net  obtained  in  practice,  and  that  only  the  solution  with  a  purely 
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hydroma gnetic  boundary  condition  on  the  elctric  field  has  physical  significance. 

It  is-  suggested  that  i  ture  research  into  tne  nature  of  hydromagnetic  shocks 
propagating  into  "cold"  upstream  states  might  profitably  include  the  effects  of 
radiative  nonequilibrium  and  photo -ionizing  reactions  in  the  analytical  models. 
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APPENDIX  -  INTEGRATION  OF  THE  NAVIER -STOKES 

SHOCK  STRUCTURE  EQUATIONS  IN  UN-IONIZED  ARGON 

The  numerical  solution  of  the  imbeii^ed  perfect-gas  viscous  shock  structure, 
while  not  entirely  straightforward,  is  well-understood  and  will  be  discussed  briefly 
here.  The  applicable  differential  equations  in  u^  and  T  are  [cf.  Eqs,  (67)  and 
(68)]; 


du  , 

—  =  Re  •  F  (u  T)  Re  [  u 

A  A  -X 

dx 


(Al) 


=  Re  •  G(u  ,  D 
dx  x 


_16 

27 


Mj* 

Re  [ T  -  1  +  —  (  u  3 
3  x 


1)  ]  • 


(A  2) 


Ordinary,  numerical  integration  of  such  differential  equations  as  an 
initial-value  problem  would  be  indicated.  It  is  well-known, however,  that  this 
is  not  possible  for  this  particular  system  because  the  derivative  dT  /  du^  = 

G  (u^  ,  T  )  /  F(ux,T)  becomes  indeterminate,  of  the  form  0 /0,  at  the  upstream 
and  downstream  states  where  F  =  G  =  0;  moreover,  the  singular-point  is  of 
the  node  type  upstream  and  of  the  saddle-point  type  downstream; 


4* 


I  iiwuWJ 
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consequently,  a  stable  numerical  solution  is  obtained  by  integrating  from  the 
downstream  state  toward  the  upstream,  but  not  vice  versa.  It  might  be 
mentioned  that  the  energy  equation  has  an  exact  integral  for  Pr  =  3/4,  in 
which  case  singular  points  in  (u  ,  T)  space  are  irrelevant  to  numerical 
integration  .  Remember,  however, that  in  the  present  problem  Pr  =  2/3  so  inte¬ 
gration  must  proceed  backward  from  the  vicinity  of  the  downstream  point. 

We  can  obtain  consistent  initial  values  for  (u^  ,  T)  in  the  neighborhood 

>Jc  # 

of  (u^  ,  T  )  ,  where  the  asterisk  (*)  denotes  the  downstream  state  of  the 
perfect  gas  shock,  provided  we  know  the  value  of  the  derivative  asymptotically 
downstream,  viz.  lim  x  — )  co  dT  /  du^  .  To  this  end,  consider  the  situation 
when  u^  and  T  are  perturbed  slightly  an  amount  A  and  £T  from  their 
downstream  values 


ux  =  Ux  +  -ux  1  T  =  T*  +  aT  , 


$  $ 

where  }ux/ux  *  <  1,  iT/T  <<  1.  To  s  implify  the  algebra,  we  make 
the  reasonable  (for  the  present  problem)  assumption  that  M*  >  >  1  so 
that  all  terms  of  order  1  /M '  or  less  compared  to  unity  will  henceforth  be 

jjt 

dropped.  The  downstream  values,  from  Eq.  (69b),  become  u^  =  1/4  , 

T*  p 

T  =5Ml/l6  and  the  corresponding  near-downptream  velocity  and  temperature 


U  =  -  +  Au 

x  4  x 


— — -  M  +  AT 
16  x 


,  'y  * 


Substituting  Eq.  (A3)  into  Eqs.  (Al)  and(A2)  and  dropping  perturbation 
terms  consistent  with  4ux  1/4  and  *T  «  5M*  /l6  yields  the 

linearized  equations 


(A  4) 


Dividing  (A  4)  into  (A5)  gives 


Now,  making  use  of  the  identities 

dT 

lim  x  - — f  oo  - 

du 

x 

and  evaluating  Eq.  (A 6)  at  x  -  oo  in  the  above  leads  to  an  equation  for  (dT/du)*  which  is 
exact  at_  the  dowmtream  sii^ular-poirt  (accepting ,  of  course,  the  approximations  related 

to  i/Mx2<<  1): 


^  . 
1UX  \i*xJ 
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and  which  has  the  two  roots 


of  which  only  the  latter  has  physical  significance. 

In  order  to  begin  numerical  integration  it  was  first  assumed,  quite 
arbitrarily,  that  we  were  at  the  point  where  A  =  1  /100  <  <  l  /4  ;  the 
consistent  value  of  the  temperature  perturbation  is,  from  the  linearized  analysis, 
A  T  ~  (dT/dux  )*•  Au^s -TM.*/  3600  .  The  initial  values  of  ux  and  T  used 
to  start  numerical  integration  at  this  point  follow  immediately  from  Eq.  (A  3); 
furthermore,  it  was  assumed  that  x  =  0  here,  in  order  to  match  the  relaxation 
zone  solution  which  (as  discussed  in  Sec.  5)  proceeds  by  forward  integration 
from  x  =  0  toward  x  — ^  +  oo  .  Integration  of  Eqs .  (Al)  and  (A  Z)  was 
carried  out  in  physical  space  by  conventional  numerical  techniques  from  x  =  0 
coward  x  =  -  oo  until  the  velocity  and  temperature  came  arbitrarily  close 
to  their  upstream  values  ux  =  T  =  1.  In  these  calculations  the  temperature- 
dspend?nce  of  the  Reynolds  number  was  given  by  R  (  T  )  =  1 . 65  M  T"3/4 
as  indicated,  for  un-ionized  argon,  by  Eq.  (60b). 
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Fig.  1 .  Sketch  of  gas -ionising  shock  structure  gepmetr>r  in  an 

oblique  magnetic  field.  The  cartesian  (x  ,y',s')  coord¬ 
inate  system  is  shock- fixed.  The  "primes"  denote 
physical  (dimensional)  quantities. 
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Fig.  5.  Nonequilibrium  structure  of  a  slow  (Mj-ZO,  gas-ionizing  shock  front, 

of  the  45°  upstream  magnetic  field  femily  (Bx  =BZ1  =  1),  propagating  into  un-ionized 
Argon  at  a  pressure  pi  =  1.0  mm  Hg  and  temperature  T' =300°K  computed  with  the 

ZND  approximation  for  Ey  =i.o.  This  case  corresponds  to  the  pure  hydromagnetic 

boundary  condition* on  the  electric  field.  The  scale  has  been  stretched  by  a  actor 

of  ten  for  .x<0  compared  to  x>0  to  show  the  imbedded  Navier- Stokes  viscous  shock  more  clearly 
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Fig.  7.  Variation  of  nondimenslonalized  downstream  streamwise  velocity 
ux  ,  transverse  magnetic  field,  BZa,  degree  of  ionization  a2  and 
nondimenslonalized  relaxation  length  tr,  for  various  shock- frame 
(and  corresponding  lab-frame)  electric  fields  for  a  slow  shock  with 
M,  =20,  Mi,  =  1  /V*,  p,  =  1 . 0  mm  Hg  and  temperature  T1  =  300°K.  A 
unique  value  of  the  electric  field  is  not  defined  by  the  structure. 
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IS.  ABSTRACT 


This  study  deals  analytically  with  th~  structure  of  gas -ionising  hydromagnetic 
shock  waves.  Since  these  waves,  by  definition,  must  have  non-electrically-con- 
ducting  upstream  states,  their  existence  at  very  high  shock  temperatures  must  be 
ruled  out  on  the  physical  grounds  that  forward- radiated  precursor  ionisation  make 
the  unshocked  gas  conducting.  A  "low  temperature"  collisicnally-ionising  shock 
with  oblique  magnetic  field  is  studied  here  to  determine  whether  certain  concepts 
which  exist  in  the  current  literature  are  relevant.  Nondimensionalised  equations 
governing  the  nonequilibrium  structure  of  such  a  front  propagating  into  un-ionized 
argon  are  formulated  using  ionisation  rates  and  an  electron  energy  equation 
developed  in  an  earlier  paper.  Comparison  of  the  magnitudes  of  viscous  and 
magnetic  Reynolds  numbers  within  this  front  indicates  that,  if  a  structure  exiSv.  , 
it  must  consist  of  a  narrow  "imbedded"  viscous  shock  standing  upstream  of  a 
much  wider  hydromagnetic  interaction  and  ionization  relaxation  zone.  Hence,  a 
modified  form  of  the  Zeldovich-von  Neumann-Doring  (ZND)  approximation  is 
applicable  to  the  structure  problem.  It  is  thown  that  iti  this  approximation  non¬ 
trivial  steady-state  structures  cannot  be  constructed  for  "fast"  gas-ionizing 
shocks.  On  the  other  hand,  solutions  are  ppssible  for  "slow"  waves,  and  these  an 
obtained  numerically  for  a  family  of  hydromagnetically  oblique  shocks  at  Mach 
number  Mi  =  20  and  Alfv4n  number  M^sl/v/7  with  parametrically  varied  values  of 
the  upstream  electric  field.  In  contrast  to  previous  expectations,  the  upstream 
electric  field  is  not  uniquely  defined  by  the  structure.  | 
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